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Abstract 



We present data of Monte Carlo simulations for monodisperse linear polymer 
chains in dense melts with degrees of polymerization between TV = 16 and 
N = 512. The aim of this study is to investigate the crossover from Rouse- 
like dynamics for short chains to reptation-like dynamics for long chains. To 
address this problem we calculate a variety of different quantities: standard 
mean-square displacements of inner monomers and of the chain's center of 
mass, the recently proposed cubic invariant [U. Ebert et al., Phys. Rev. Lett. 
78, 1592 (1997).], persistence of bond-vector orientation with time, and the 
auto-correlation functions of the bond vector, the end-to-end vector and the 
Rouse modes. This analysis reveals that the crossover from non- to entangled 
dynamics is very protracted. Only the largest chain length N = 512, which 
is about 13 times larger than the entanglement length, shows evidence for 
reptation. 
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I. INTRODUCTION 



A detailed description of the slow dynamics of entangled polymer chains would be of 
great interest for better understanding the physical properties of dense polymer melts or 
solutions. Several attempts to solve this problem have been made. In a seminal paper de 
Gennes introduced the concept of reptation [|TJ. The fundamental idea is that the motion 
of a polymer is not spatially isotropic, but has to occur along the contour of a "tube" 
which is formed by the surrounding chains. These surrounding chains are treated as a 
fixed network of impenetrable obstacles. The simultaneous motion of the polymers and 
possible correlations resulting from that are disregarded. The elaboration of this idea has 
led to a rich scenario of testable predictions p|-§J . Critical analysis of these predictions by 
experiments and computer simulations fT(l|-|T2]| suggest that reptation is an important 
relaxation mechanism. However, it cannot provide a quantitative explanation of all dynamic 
features observed. These deviations could result from additional relaxation processes which 
blur the clear signature of reptation for typical experimental chain lengths N and become 
negligible only in the asymptotic limit of extremely large N (see, however, [lj] for a different 
asymptotic behavior). 

These additional relaxation mechanisms could involve fluctuations of the tube's contour 
|I^TB|, reorganization of the tube |HJ by release and creation of topological constraints 
|5| JT^JT8[| , elastic distortions of the entanglement network |19[ or long-lived density fluctu- 



ations pil. All these extensions assume the existence of a tube. None of them explains 



its microscopic origin. An attempt to develop a microscopic theory was made by the de- 
velopment of a mode-coupling theory for polymers [pTH 26|] . This approach leads to similar 
predictions as reptation theory without invoking a tube hypothesis. 



However, some doubt was cast on this concept by simulations of a model by Shaffer [27 



which respects excluded volume interactions between the monomers of the chains, but al- 
lows Monte Carlo moves that lead to intersections between bonds connecting the monomers. 
For long chains in this model a simple Rouse-like behavior is found. Only when the condi- 
tions that bonds must not intersect is used as an additional constraint in the dynamics, a 
reptation-like behavior results from the simulation. These results, where static correlations 
are precisely identical for both models, but one dynamic version exhibits Rouse behavior and 
the other reptation-like behavior, seem to imply that the topological constraint of the non- 
crossability of chains is not an automatic consequence of the excluded volume interaction. 
In this context, we note that many other attempts have been made to account for the slow 
dynamics of entangled polymer chains by concepts differing from reptation, e.g., collective 
motion of many entangled chains ||28|| , but we are not attempting to critically evaluate all 
these theories here. 

In view of these diverse attempts to improve on reptation theory, it is also beneficial to 
consider a model which exhibits pure reptation behavior in the asymptotic limit of large 
chains and long times, and for which finite- N and finite-time corrections can be calculated. 
Evans and Edwards proposed such a model in which a single chain moves through a reg- 
ular array of impenetrable obstacles p9|| . The obstacles do not influence the equilibrium 



configuration of the chain, but only the accessible relaxation moves of its monomers. This 
model has recently been reanalyzed analytically p3| . |3T|| . Detailed comparisons with com- 
puter simulations [K0,^,Q reveal an extremely slow crossover to the asymptotic power laws 
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characteristic of reptation. Even in the ideal situation, where the tube diameter is equal to 
the bond length, the power laws only become clearly pronounced for iV > 100. If the lattice 
of obstacles is diluted, the tube diameter increases and the crossover to reptation is shifted 



to larger chain lengths |33 



The latter study resembles more closely the situation encountered in computer simula- 
tions of melts where the tube diameter is in general larger than the bond length. Taking 
into account that simulations typically work with iV < 1000, the analysis of the Evans- 
Edwards model suggests that the results of polymer melt simulations are characteristic of 
the crossover from non-entangled, Rouse-like dynamics for short chains to (slightly) entan- 
gled, reptation- like dynamics for long chains. With the present study we want to further 
investigate this crossover by a Monte Carlo simulation of the bond-fluctuation model. Our 
paper is organized as follows: Section [TT] introduces the model and describes some of its static 
properties. Section |TXT| represents the main part of the paper. It deals with dynamics of the 
model. Here, we compare various mean-square displacements and orientational correlation 
functions, extract the chain length dependence of relaxation times and of the diffusion co- 
efficient, and estimate the parameters of reptation theory, such as the tube diameter or the 
entanglement length. The final section |V| summarizes our conclusions. 



II. MODEL AND STATIC PROPERTIES 

The aim of the present work is to extend and to complement previous studies of the 
crossover from Rouse to reptation dynamics using the bond-fluctuation model H3"4H36|. The 
bond-fluctuation model is a coarse-grained lattice model |T^J3^rj3|, in which a polymer is 



represented by a chain consisting of N monomers. Each monomer occupies a unit cell of a 
simple cubic lattice (unit cell = 8 lattice sites). This increased monomer size allows the bond 
vectors I to vary in both length and direction to a much larger extent than in simpler lattice 



models where a monomer is associated with a single site |37j. In total, there are 108 bond 
vectors [5!| giving rise to 87 bond angles and 5 different bond lengths (2, a/5, a/6, 3, a/To). 
The average bond length and bond angle depend on the volume fraction, <fi = 8NP/L 3 , 
of occupied lattice sites (P = total number of polymers in the simulation box, L = linear 
dimension of the box). Choosing = 0.5 the model realizes typical properties of dense melts 



34[ . For this density the mean bond length and mean cosine of the bond angle are given by 
/ = (Z 2 )V2 = 2.64 and a = (cosfl) = -0.1055, respectively gO]. 

To analyze the change in dynamic behavior with increasing N simulations were performed 
for different chain lengths ranging from iV=16toiV = 512 (see Table |). This extends 
previous studies |34J by a factor of 2 only, but it amounts to an increase in relaxation time 



by almost an order of magnitude (cf. Table |). Furthermore, the system size could also be 
chosen much larger. For most chain lengths the linear dimension of the simulation box is 
more than a factor of 2 greater than the average end-to-end distance R e , i.e., L > 2R e (with 
the exception of N = 512 where L = 128 ~ 1.75i? e only; see Table [|). Since finite-size effects 
may be expected if L < 2R e , Fig. [I] compares the mean-square end-to-end distance and 
mean-square radius of gyration R^ of the present study with the values of earlier works 
3"4] , |3"5| . Two conclusions can be drawn from this figure. First, finite-size effects did not 



affect previous (static) results. Within the statistical uncertainties the old data for N = 200 
(R e ~ 45, L = 40 p4|) agree with the present analysis. Second, the bond-fluctuation model 



3 



approaches the asymptotic ideal chain limit, i.e., B% = 6R 2 = Cool 2 (N — 1), if N > 100. 
Here, (= hin^^oo B%/[(N — l)/ 2 ]) is the characteristic ratio which measures the swelling 
of a chain due to its intrinsic stiffness compared to a pure random walk. Using ~ 1.56 
(see Fig. |l|) and l v = 1(0^ + l)/2 pHl , the persistence length, / p , of the model is given by 
l p ~ 3.38. 

The shape of a polymer is not spherical. A good indicator of its geometrical form is the 
gyration tensor Q whose elements are given by f4"2|-pB| 

1 N 

Qaf3 = ^^(-Ra.n ~ Rcm,a)(R/3,n ~ ^cm.jj) («, /3 = 1, . . . , 3) . (1) 

n=l 

Here, R a , n and R C m,a denote the ath-component of the position vectors to monomer n and 
to the chain's center of mass, i.e., R n = (-Ri, n , R2,n, R3,n) and R cm = (R C m,i, R C m,2, -R C m,3)- 



The trace of this tensor provides another way to calculate the radius of gyration |42| 

(TrQ) = (\ 1 ) + (\ 2 ) + (\ 3 ) = R 2 g , (2) 

where X a are the eigenvalues of Q. If the shape of a polymer was spherical, one would have 
(Ai) = (A2) = (A3). Contrary to that, Fig. § shows that the biggest eigenvalue, (Ai), is much 
larger than the other two eigenvalues. This implies that the shape of a chain is distorted 
with respect to a sphere: It is stretched in direction of the long principal axis and shrunk in 
directions of the axes corresponding to (A2) and (A3). Therefore, a chain rather resembles a 
"flattened American football"[]. This characteristic property can be exploited to construct 
an efficient coarse-grained model for dense polymer melts P7LIIS . 



In the course of the simulation, a monomer and a lattice direction (out of 6) are selected 
at random, and a jump by one lattice constant is attempted in that direction. The jump is 
accepted if the resulting bond vectors belong to the allowed set of bonds and if the targeted 
lattice sites are empty. Otherwise, the move is rejected. This random displacement simulates 
the impact of a stochastic force which is exerted on a monomer by its local environment. 

In the single chain limit, this local jump model is expected to give rise to Rouse dynamics 
if hydrodynamic interactions are absent ||]. Since hydrodynamic forces are mediated by 
solvent molecules, an isolated chain of the bond-fluctuation model provides an example for 



this conjecture [p!2|,p6|. However, if the volume fraction increases, deviations from Rouse 



dynamics are observed due to the mutual interaction of the chains and due to the onset of 



entanglements [12|]. Nevertheless, the Rouse model is generally believed to provide a reliable 
description of short (non-entangled) chains in a melt. Therefore, we determined its basic 
variables, the Rouse modes 



x P (t) = ±f2 R ^ cos in b = M,...,iV-i), (3) 

n=l 



Recent studies [^5 .46] of the shape of Gaussian chains indicate that the visualization of a random- 
walk polymer as a flattened American football is not completely correct. The density distribution 
of monomers in the coordinate system of the principal axes exhibits a slight minimum at the origin 
(i.e., at the center of mass) for the largest axis, whereas it has a maximum for the other two axes. 
Therefore, the shape is rather dumbbell-like. 
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where R n {t) is the position of the nth monomer at time t. At t = the (static) correlation 
of the modes reflects the structural properties of the polymers. For ideal random walk chains 
it is given by [pj 



(X p (0) ■ X fl (0)> = 6. 



pq 



IN 



sm(pn/2N) 



p/iV<l 



1 Nb 2 



(for p > 0) . (4) 



In this equation b denotes the effective bond length defined by b 2 = C^l 2 0. However, we 
use b 2 = R 2 /(N — 1). This ratio is (slightly) chain-length dependent. It varies from about 
8.7 for N = 16 to about 10.6 for N = 512, i.e., by about 18%, in our simulation. 

Using this identification Eq. suggests that a plot of 8N(N — 1) (X p (0) 2 ) / R 2 e versus 
p/N should yield a master curve. Figure shows a test of this prediction. In fact, the data for 
all chain lengths collapse onto a common curve which is nicely described by [sin(p7r/2iV)]~ 2 
if p/N < 0.05. However, if p/N increases, deviations between the simulation results and 
the Rouse prediction gradually develop. The Rouse theory overestimates the correlation, 
especially for p/N > 0.3. These large modes probe local distances along the backbone of a 
chain. For instance, p/N > 0.3 corresponds to subunits of a trimer and smaller than that. 
On these local scales, the intrinsic stiffness of a polymer should be taken into account. 

The simplest way to achieve this consists in replacing the random walk by a freely 
rotating chain flj4l||. This is a chain model, in which the mean bond length and bond angle 
are fixed and each bond is allowed to rotate freely around the direction of the preceding 
bond. Starting from 



1 N 

(X 2 p (0)) = — (# n (0)-i* m (0))cos 



pm{n — 1/2) pn(m— 1/2) 

COS 



n,m=l 



N 



N 



this implies to write (remember a = {cos 8 

n—l m—l 



-0.1055) 

n— 1 m—l 



(R n -R m ) = j2J2^-^ = l2 J2Y,(- 

i=l j = l i=l j = l 

instead of (U ■ lj) = l 2 5ij for a random walk. Using furthermore 



a 



N 

E 



pn(n — 1/2) pn(m— 1/2) 
cos — cos - 



n,m=2 

N n N 



N 



N 



n—l m—l 
i=l j=l 



\i-3\ 



pTxin — 1/2) pixim — 1/2) 

1^1^2^ cos cos 



n=2 i=1 m=i 

one obtains 



N 



N 



-a 



iTTl— 1+1 



-a) - 1 



1 (-a 
- + ^— r- 



.n— i+1 



a) - 1 



Rl 



sin(pvr/2A^) 



4(-a) 



1 - 2(-a) cos(j9vr/A^) + {-af 



1 + 



1 2a(l + a) 



sin(p7r/A^) 



N I- a 1 + 2acos(p7r/iV) + a 2 tan(p7r/2iV) 



(5) 
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where the ~-sign indicates that \a\ = 0.1055 <C 1 was used to write the 1/iV-correction in 
the second line of Eq. (|5|) in the present form. Although this correction violates the p/N- 
scaling, its numerical value is always smaller than 0.03 (= "worst case" for p = 1 and N = 16) 
and thus negligible. Therefore, we compare the first line of Eq. (|5]) to the simulation data. 
Figure [3] shows that the approximate consideration of stiffness by Eq. ([5]) reasonably accounts 
for the suppression of correlations below the Rouse prediction Similar approaches are 
pursued in ]50 , CT and tested against experiments and simulations [ p^ , ^ ]. 

Figure |] tests another prediction of Eq. the orthogonality of the Rouse modes, i.e., 
(X p (0)X q (0)) ~ S pq , for two chain lengths N = 16 and iV = 128. Whereas the auto- 
correlation of the first mode is 1 to 2 orders of magnitude larger than cross-correlations 
with p > 1 for both chain lengths, this difference in magnitude is only preserved for the 
short chain when increasing p towards p = N — 1 (see Fig. |]b). On the other hand, cross- 
correlations of the same amplitude as the self-correlation develop between p = 127 and the 
smallest modes for N = 128. If one interpretes a difference of 1 to 2 orders of magnitude 
between self- and cross-correlations as a numerical realization of (X p (0)X q (0)) ~ S pq , the 
full spectrum of Rouse modes is only delta-correlated for N = 16, while small and large 
modes interfere for large chain lengths. 



III. DYNAMIC PROPERTIES OF THE MELT 

This section discusses the simulation results for dynamic properties of the polymer melt. 
It is split into three subsections. The first deals with an analysis of various mean-square 
displacements. The presentation of the data is very much influenced by a recent thorough 



theoretical and numerical study of the Evans-Edwards model [p0| -|33|| . The second subsection 
presents the relaxation behavior of the correlation function of the bond- vector, of the end-to- 
end distance and of the Rouse modes, whereas the last subsection discusses the chain-length 
dependence of the corresponding relaxation times and of the diffusion coefficient. 



A. Mean-Square Displacements 



Important quantities to study the crossover from Rouse to reptation dynamics are the 
following mean-square displacements: 



9i(t) 
92(t) 

9s(t) 
9*{t) 



-Rat/2 (*) — -Rat/2 (0) 

■Rat/2 (t) — -Rat/2 (0) — R cm (t) + -R cm (0) 
1 2 \ 

Remit) — -R cm (0) ) , 

Rt{t) - R 1 (0)] 2 ) + ( \R N (t) - R N {0) 



(6) 



Here, gi is the mean-square displacement of the inner monomer of a chain (situated at 
position R^/ 2 it) at time t), g 2 is the same displacement measured relative to the motion 
of the chain's center of mass (-R cm (t) is position of the center of mass at time t), g 3 it) is 
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the mean-square displacement of the center of mass, and g± is that of the end monomers. 
Within the framework of reptation theory one expects |pj^,P,12 



9i(t) 




if t < r e , 

if r e < t < r R , 

if r R < f < r d , 

if r d < i , 



03 (*) 



t/JV 

(t/N 2 ) 1 / 2 
t/N 2 



if t < r c , 

if r e < i < r R , 

if r R < * . 



(7) 



At early times, a chain does not feel the entanglement constraints imposed by its neighbors. 
Assuming the Rouse model to provide a realistic description for non-entangled chains, one 
expects gi ~ t 1 ^ 2 and g^ ~ t in this case. If t equals the entanglement time r e , the constraints 
begin to dominate the polymer dynamics. For larger times, a chain moves as if it were 
confined in a tube created by its neighbors: It can only slide along the tube axis (the so- 
called "primitive path"), whereas motion perpendicular to it is suppressed. Since the tube 
can be thought of as an envelope around the chain's random-walk-like configuration, the 
chain performs Rouse dynamics along a random walk. This leads to g\ ~ t 1 / 4 and g% ~ t 1 / 2 . 
If t reaches the Rouse time tr, the chain configuration is relaxed inside the tube. For larger 
times, the center of mass diffuses freely and the inner monomer begins to diffuse out of the 
tube. The monomer leaves the tube completely if t equals the disentanglement time r<j. 

The power laws of Eq. (|7|) represent the asymptotic behavior. They can be observed 
clearly only for very large N, where the time scales r e , tr and r d are well separated. For 
smaller N significant corrections must be expected. This has recently been demonstrated by 
a detailed theoretical analysis |31| of the Evans-Edwards model p9| and by comparing the 
analytical results with the outcome of simulations [30,32|. The Evans-Edwards model is a 
lattice model, in which a single chain moves on a simple cubic lattice through impenetrable 
obstacles placed in the centers of each unit cell of the lattice (a monomer occupies one 
lattice site). These obstacles impose a strong geometric confinement on the motion of the 
chain: The tube diameter di is of the order of the lattice constant. Even under these 
very favorable conditions long chains and long times are required to display the asymptotic 
power laws clearly p0|j3^| . If one wants to analyze the motion of the inner monomer in a 
different model, the experience with the studies of f30|-|33| suggests to use the following three 
quantities: the mean-square displacement #2 which behaves as 



92® 



gi(t) if t < r d , 



R 2 if r d < t 



(8) 



so that a transition from ~ t 1 ^ to ~ t 1 ! 2 cannot be interpreted as a crossover to free 
diffusion due to the saturation of #2 at late times, the "cubic invariant" [|30| -33 



i 



[ R <*,N/2(t) - R a , N /2(0) 



a=l 



9i(t) 



(9) 



which is predicted to be much less plagued than gi by preasymptotic corrections (for the 
Evans-Edwards model), and the ratio g&(t) / gi{t). This ratio should start around 1 at early 
times, develop a maximum at intermediate times and approach 1 if t — *■ oo (= free diffusive 
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limit for both g± and g\). For Rouse dynamics the maximum occurs at t < tr and has the 
value g^/gi = 2. If the chain reptates, the analysis of the Evans-Edwards model suggests 
that the position of the maximum is still given by t ~ tr, but its amplitude should be much 
larger than 2 (g 4 / gi w 4^ for iV -► oo) U^H. 

Figures ^]and|^ depict the time dependence of and #2j respectively. The ordinates are 
divided by t 1 ^ to highlight the expected asymptote for the onset of reptation. The figures 
show that the initial increase of g§ and g 2 , which is roughly compatible with a t 1 / 2 -behavior 
for 200 < t < 5000, considerably slows down when t > 10 5 . However, an indication of 
the t 1//4 -power law is only visible for N = 512. Even for this chain length there is no clear 
sign of a subsequent t 1//2 -increase, but only a fairly protracted transition to the large-time 
limits g& ~ t and g 2 — > i? 2 . All other chain lengths already cross over to these limits for 
t > 10 5 (this is especially well visible for g 2 ). Therefore, much larger (or stiffer P5,pl,p3|) 



chains and still better statistics are needed to clearly separate the different power laws and to 
eventually distinguish between reptation theory and alternative approaches, such as polymer 
mode-coupling theory fl21-2"6|. 



Nonetheless, estimates of the tube diameter g?t and the entanglement time r e can be 
obtained from the figures by posing 

d 2 

9q{t c ) ^ g-i. (r e ) ~ gi(r e ) = y , (10) 

where r e is defined by the intersection point of the t 1 ^ 2 - and t 1//4 -power laws. The factor 1/3 in 
Eq. ([T0| ) was proposed in |56| , |5"7| and justified by the following argument: The tube diameter 



can be interpreted as the end-to-end distance of a chain segment with iV e monomers. This 
is the largest segment which does not feel the entanglement constraint yet. It relaxes as 
a Rouse chain. So, r e = tr(A^) which gives gi(r c ) pa <i 2 /3 |§]. Qualitatively, the factor 
1/3 also appears reasonable. If we assume that the end-to-end vector d 2 ^ is predominantly 
oriented along the tube axis, the average extension of the tube perpendicular to its axis 
should be proportional to the radius of gyration d\/Q and thus smaller than d\ ||33[| . 

If one accepts these arguments, Eq. ( |T0"D yields d T pa 17.9, 19.3 and r e w 84053, 100175 
for g 2 and g^, respectively, so that on average d^ ~ 18.6 and r c pa 92114. Using the result 
for g?t we can furthermore estimate the entanglement length iV e by requiring d\ = R 2 (N e ). 
To this end, we fitted the simulation data for R 2 in the range 16 < iV < 64 by a power law, 
which gives R 2 {N) = 6.872iV L08 so that iV c pa 38. This value is about 26% larger than the 
original estimate f3~4}] , but coincides very well with the result of a recent study 

Figure |^ shows the simulation results for the ratio of the mean-square displacement 
of the end monomer and the inner monomer, gi{t) / g\{t). The ratio starts around 1.5 at 
early times, exhibits a maximum at intermediate times and tends to 1 for large t. For 
iV < 64 the maximum occurs close to (74/(71 = 2, i.e., to the value expected from the 
Rouse model, whereas it is larger for iV = 128 {g^/gi ~ 3). An increase of the maximum 
beyond the Rouse limit with growing chain length is predicted theoretically for the Evans- 
Edwards model and proposed as a sensible indicator of reptation dynamics [0-33]. For 



a reptating chain the disparity between g± and g\ has to increase because the motion of 
the inner monomer is strongly confined by the tube at intermediate times, whereas end 
monomers always take part in tube renewal. Our data for iV = 128 are indicative of such 
an enhancement beyond the Rouse limit and are qualitatively similar to simulation results 
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of the Evans-Edwards model with a randomly dilute obstacle lattice (see Fig. 5 of [|33|| ). 
Unfortunately, we have to restrict the discussion to iV < 128 because the calculation of the 
ratio makes the statistical inaccuracies become so pronounced that a quantitative analysis 
for N = 512 was not possible. 

In addition to the displacement of inner and end monomers the motion of the center of 
mass is also interesting. Figure |] plots (?3(t)/t 3//4 versus time for all chain lengths studied. 
This representation of the ordinate was motivated by Eq. (|?J). If a chain reptates, g3(t)/t 3 ^ 
should exhibit a minimum close to t — tr, which is preceded by a t _1//4 -descent and followed 
by a t 1//4 -ascent. In analogy to #2 and g6, a signature of this behavior emerges for iV = 512 
only, whereas the transition to free diffusion intervenes for smaller chain lengths. 

The data reveal that gs is subdiffusive at early times. Particularly for iV = 512, where r e 
and the crossover-time to free diffusion are well separated, this behavior is very pronounced. 
The mean-square displacement increases as g% ~ t x with an effective exponent x ~ 0.84. 
Subdiffusive motion of the center of mass at early times is not unusual. It is generally found 
in computer simulations of linear polymers (see for reviews and [58] for a chemically 

realistic model of polyethylene) and of rings ||59|| , and it has also been observed lately in 
experiments (on short chains) ||53|| . The subdiffusive behavior should be considered as a 
condensed-phase effect. In a dense melt, a polymer is intermingled with many other chains. 
Their presence impedes the motion of the center of mass of the tagged chain, leading to 
subdiffusive behavior. Evidence for this interpretation is provided, for instance, by simu- 
lations at different densities. If the density decreases towards the single-chain limit, the 
bond-fluctuation model yields g% ~ t for all times (i.e., typical Rouse behavior) ||T2| , |34] , |36 



Another evidence stems from recent molecular dynamics simulations of the Kremer-Grest 
model [3C]. In this study, the velocity-autocorrelation function of the center of mass was de- 



termined. The correlation function becomes negative at intermediate times and approaches 
from below. This means that the surrounding polymers exert a force on the tagged chain, 
which reverses its velocity and tries to push it back to its original position. 

Since the subdiffusive motion of g^ is already present for short non-entangled chains, 
it should replace the expectation from the Rouse model g 3 ~ t for t < r c [see Eq. (0)]. 
Therefore, one can try to determine r e from the intersection of g% ~ t a84 and g% ~ t 1 ^ 2 . This 
yields r c w 268718, which is about a factor of 2.5 larger than the results derived from g e or 
gi. Nevertheless, if one accepts this value, another estimate of can be obtained by posing 



b 2 d\ 



3R 2 e \(b 2 



1/2 



(11) 



where the factor 1/3 stems from Eq. ( |10|) and the requirement gi(r d ) = ^(ra). In Eq. (pT|), 
( denotes the monomeric friction coefficient and r e is (assumed to be) given by [|1^ 



so that 



dn 



d'Y 

T 



3#3 



knT 



)Ri 



,1/4 



(12) 



(13) 



Using g 3 (r c ) w 10.37 and R 2 = 5348, Eq. (TTJ) leads to d T w 20.3. This value is compatible 
with previous estimates from g 6 and g^. However, due the small exponent 1/4 in Eq. ([13]) 
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the result for the tube diameter is not very sensible to variations of r e . For instance, when 
using the average r e from g 6 and g 2 , i.e., r e 92114, so that gy,ij e ) ~ 4.27, one obtains 
d<i 16.2 which is also in reasonable agreement with the results derived from the monomer 
displacements. 

On the other hand, the power law g% ~ t 1 ^ 2 provides a further possibility to determine 
c?t by the intersection with the diffusive behavior g^ ~ t. The intersection point should 
define the Rouse time tr (« 4.34 x 10 6 ), and one expects the value ^(tr) to be close to 
d\. If we determine the prefactor by the same argument which led to Eq. (|TT|), we obtain 
93{ t r) — d\/?>. Using ^(tr) « 54.8, we find d T « 13. This value seems too small in 
comparison with the previous results which indicates that the Rouse time is presumably 
larger than determined. Therefore, tr = 4.34 x 10 6 rather corresponds to a lower bound 
than to a reliable estimate. 



B. Correlation Functions 

If a chain starts to reptate, the motion of the monomers perpendicular to the primitive 
path is restricted to displacements smaller than the tube diameter. This implies that the 
decay of orientational correlation functions of chain segments should significantly slow down 
if the size of the segment becomes comparable to dy. Therefore, correlation functions which 
probe reorientations on different length scales along the backbone of a chain are interesting 
quantities. 

The smallest segment along the backbone of a chain is the bond vector. The correlation 
function of the bond vector can be defined by 

1 JV ~ 1 

Mt) = {N _ 1)P E (kW • z «(°)) > ( 14 ) 



n=l 



where I 2 = (Z 2 (0)) is the mean-square bond length. This quantity is also interesting because 
it should be closely related to the shear modulus f26|. Figure || illustrates the time- and 
chain-length dependences of 4>b{t). The data are presented in a scaling plot where the time 
axis is given in units of the relaxation time, Tb, of the bond vector defined by 0b(i"b) = 1/e. 
This time scale depends weakly on chain length (see Table |). If Tb was the only relevant time 
scale, all 0b(t)'s, measured for different N, should collapse onto a master curve. Figure |9| 
shows that such a collapse is realized for about 75% of the relaxation, whereas the final 
decay to zero depends on chain length. Larger chains relax more slowly than shorter ones. 

Such a behavior is already expected from the (discrete ||49|| ) Rouse model which expresses 
(j) h (t) as 



1 N ~ 1 1 JV ~ 1 t 2 

Mt) = ^tE Mt) = j^j £ exp ( - A) , (15) 
p=i p=i K 

where tr is the Rouse time and <f> p {t) the correlation function of the pth Rouse mode @J49[ 

(X p (t) ■ X p (0)) ( tp 2 x 

Mt) = {xm) = exp ( - —) . (16) 
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If t <C tr, the relaxation behavior of the sum in Eq. (|I~5"D is dominated by large p's, which 
yields 0b (t) ~ l/t 1//2 independent of A" (in the continuum limit). On the other hand, if 
t > tr, 0b (t) relaxes the more slowly, the longer the chain length. These predictions seem 
to be in accord with the simulation results. However, there are differences at short times 
for all chain lengths and at long times for large N. 

At short times, the Rouse formula underestimates the relaxation strength of 0b (t). 
This discrepancy is caused by Eq. (|TB]) which assumes an exponential decay for all Rouse 
modes. However, especially the simulation results for large Rouse modes do not decay 
in a simple exponential fashion (see discussion of Fig. |T2|). These large modes determine 
the behavior of 0b (t) at short times, whereas the late-time relaxation is dominated by the 
smallest Rouse mode which is almost exponential for A" < 128. Nevertheless, Eqs. (|i~5"D 
and ( |T6| ) only provide a reasonable description for the A^-dependent tail of <fih(t) if N = 16. 
For larger A" the simulation data relax more slowly and finally develop a power-law decay 
if N > 128. The power law starts around t ps r e and extends over about 3 decades up to 
the end of our simulation. With the present simulation results we cannot decide whether 
the power law is an indication of a possible two-step relaxation which would become clearly 
visible for longer chains, as suggested by polymer mode-coupling theory f24] , |2"6"|| . 



Qualitatively, the slow power-law decay of 0b if) implies that some bonds lose the memory 
of their original orientation only very gradually. Persistence of bond orientation is also an 
important notion in the development of the continuum theory of reptation |4j]. There, a key 
quantity is the probability ip{s, t) that a segment s of the primitive chain is still in the tube 
at time t. This probability is related to the projection of the (unit) tangent vectors, u(s', t), 
of all segments s' onto one primitive path segment s, i.e., onto u(s, 0) (see Eq. (6.45) of 
Ref. 0). In analogy to i/;(s,t), we define 



En>=l(ln>(t) ■ l n (0)) 



which represents the probability that the original orientation of the bond vector l n (0) is still 
present at time t. This quantity can be calculated by the (discrete fH| ) Rouse model. The 
result is 



, U s E P =i(odd)[ sin (^ 7r / iV )/ tan (W2A^)]exp(-p 2 t/r R ) 

0nW = ^v^i — ( /m/ , f , om ( 18 ) 

E P =i (odd) sin(np7r/iV)/ tan(p7r/2iV) 



p/N<l 4 



4 1 

- ^ - sin(np7r/A0 exp ( -p 2 t/r R ) . (19) 

(p odd) 

The continuum limit (|i~9|) is identical to the reptation formula for ip(s, t) if tr is replaced by 
the disentanglement time r<j. Numerically, the discrete and continuum results for n (i) are 
almost indistinguishable. Therefore, we compare Eq. (|18D with simulation data for A^ = 16 
and N = 128 in Fig. [10]. 

If N = 16, Eq. ([Tj^) provides a fairly accurate description of the simulated n (t). Both 
theory and simulation yield bell-shaped curves whose amplitudes become smaller as time 
increases. This behavior can be rationalized as follows: Since chain ends are more mobile 
than inner monomers, bond vectors close to the ends decorrelate more quickly than those in 
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the inner part of the chain. The decorrelation thus propagates from the chain ends towards 
the middle monomer. The bond vector orientation of the inner monomer is rather long- 
lived. Even at t — r ee , where the end-to-end vector correlation has decayed to about 30% 
of its original value (see Eq. (^l|) and Fig. |ll]), 4> n {t) is still about 0.5 for n = 8. This slow 



relaxation is already present for short, Rouse-like chains and should not be considered as 
a characteristic feature of reptation-like dynamics, since the Rouse and reptation formulas 
agree with one another [Eq. (|I5|)]. 

However, the simulation results for <p n {t) exhibit differences between short and long 
chains. As Fig. |H] shows, the relaxation of inner monomers becomes strongly retarded 
compared to the theoretical (reptation) prediction il N — 128. More than 50% of the 
initial orientation of the inner third of the chain is preserved even at r ee . These long-lived 
correlations, which are not present for small N, should be responsible for the power-law 
behavior of <fib(t) (see Fig. |D because they begin to develop in the same time window as the 
power law, i.e, for times larger than the entanglement time (r c /r cc « 0.04 for iV = 128), and 
last until t > r ee (the time t/r h = 10 4 in Fig. |9] corresponds to about t/r ee ~ 2 in Fig. [Ujb). 

Furthermore, they should also become visible in the motion of the monomers, particularly 
when comparing the mean-square displacement of the end monomer, (t), with that of the 
inner monomer gi(t). Figure suggests that the ratio g±{t)/gi(t) for t « r ce is much larger 
for iV = 128 than for N = 16. Our simulation data for g^/gi support this expectation (see 
Fig.©. 

If the chain length increases, one may expect the slowly relaxing zone around the middle 
monomer to grow and to influence the dynamic behavior of the chain more and more. For 
instance, one might speculate that the end-to-end vector relaxes in two steps: There is a 
first step initiated by the chain ends. Their fast motion decorrelates adjacent bond vectors 
very efficiently. However, this decorrelation does not propagate homogeneously along the 
polymer backbone towards the inner monomer as it is the case for short chains, but slows 
down in the inner part of the chain. The relaxation of this part is responsible for the second 
step. 



Some evidence for this conjecture is given in Fig. [IT] which shows the time-dependence of 
the end-to-end vector correlation function, <f) e (t), for all chain lengths studied and compares 
the simulation data with the Rouse (or reptation) prediction 

^ = { ^m )) = t^^ t,T " h (20) 



(p odd) 



where the relaxation time r ee is defined by 



(ree) = V — exp(-p 2 ) = 0.298221 . (21) 

/ — ' p Z 7l Z 
P=l * 



(p odd) 



Note that Eq. fl2"U| ) is the same for both Rouse and reptation theory, the only difference 
being that r ee is the Rouse time, tr, in the first case and the disentanglement time, Td, in 
the latter case, i.e., 
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(22) 



Td 



^&k N3 if N>>N ° 



where prefactors 1/37T 2 and 1/tt 2 are taken from M (they have been set equal to 1 in the 
discussion of the mean-square displacements in Sect. [Ill A|) . 

Equation ( |2"0"D suggests that simulation data for (j) e {t) depend on N only via r ee . By using 
a rescaled time axis, t/r ee , it should be possible to superimpose all data onto a master curve 
which is given by Eq. (|20|) . Figure [TT] confirms this expectation if N < 128. For N = 512, 
however, the relaxation is more complicated. It seems to occur in two steps. The first step 
starts at about t w r e , and the second takes place in the time window where crosses 
over to free diffusion, i.e., on the scale of the Rouse time r R (see Fig. ||). This behavior is 
similar to that expected from polymer mode-coupling theory [f24| -|26|l . The theory predicts 
two power laws: 1 — (j) c {t) ~ t 9 / 32 for r c C f < r R and 1 — e (t) ~ t 3 ^ 8 for r R t <C r RR 
(trr — renormalized Rouse time < ra). Our data are compatible with a t 9//32 -behavior for 
about one decade (0.002 < t/r cc < 0.1) and eventually with t 3//§ close to r ee . Of course, 
this can only be considered as an indication that these predictions might be relevant for 
our model. Clarification of this point requires a much better separation of the relevant time 
scales and thus simulations of longer chains. On the other hand, it is not obvious that 
Eqs. (^) and (^l|) are accurate for chain lengths in the crossover regime from Rouse to 
reptation behavior, although they describe both limiting cases of the simple Rouse and the 
"pure" (i.e., asymptotic) reptation theory. Certainly, a quantitatively reliable theoretical 
prediction for <p e {t) as well as for r ee in the regime where N and N e are comparable would 
be very desirable. 

The observed two-step-like relaxation behavior is not limited to <p c (t). It can also be 
observed for the correlation function of the Rouse modes <p p {t). Figure [12] shows scaling 



plots of p (t) for N = 1Q and N = 512. As suggested by Eq. (|16|), the scaling time, 
t p , is defined by P (r p ) = 1/e. The figure illustrates that the first Rouse mode is almost 
exponential for N = 16, whereas the relaxation of higher modes becomes progressively non- 
exponential with increasing p (smaller distance along the chain backbone). However, the 
curves cannot be described (completely) by a stretched exponential. This would imply that 
the fast decay of, say, p — 15 for <p e {t) > 1/e entails a correspondingly slower relaxation 
if (f) e (t) < 1/e. Since the data splay out with increasing p at short times, they should also 
splay out at late times, but in reverse order. There is no indication of that behavior in 
Fig. [L^a. The interplay of chain stiffness and local excluded volume forces gives rise to a 
more complicated than stretched-exponential time dependence. The same behavior is also 
observed for a bead-spring model of short chains [^TJ, and, to some extent, for iV = 512 in 



Fig. |12ib. For large chains, however, there is another factor which may contribute: cross- 
correlations between Rouse modes. Figure |3| revealed that these can be as pronounced as 
the self-correlation if p = N — 1. Since <f> p=1 (t) exhibits again the signature of a two-step 
relaxation, such as (fi e (t), this behavior could also influence the decay of p =5n(t). 
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C. Correlation Times and Diffusion Coefficient 



The discrete Rouse model yields the following expression for the relaxation time r p [JE 

7 n 2 



7p 



12 fc B T 
2 A^C /v2 



sin(p7r/2iV) 



(23) 



3 k B T 



(Km , (24) 



where Eq. (|J) was used to obtain the last line. The unit of t p is set by the monomeric friction 
coefficient ( which should be independent of the mode index p and chain length N. Within 
the framework of the Rouse model ( can be determined from the diffusion coefficient, D, of 
a chain by 

ND = ^ . (25) 

If one furthermore identifies the effective bond length with b 2 = R 2 / (N — 1) (as already done 
in Fig. IS]), Eq. (|23| ) suggests to construct a master curve by plotting 12DN(N — l)r p /R 2 
versus p/N. Figure [13] shows that such a scaling is only possible if N < 64. For these chain 
lengths the scaling yields a curve which is independent of N and coincides with the Rouse 
prediction for p/N < 0.06. On the other hand, there are deviations for p/N > 0.3, where 
Eq. ( P3"D overestimates the simulation results. 

Qualitatively, this behavior is identical to that observed for (X 2 (0)) in Fig. |3|. This 
suggests that the discrepancy between simulation and Eq. (|23|) might be related to local 
stiffness effects and that it could be removed by expressing (Xp(0)) via Eq. (|5|) or by inserting 



the simulation results for (Xp(0)) in Eq. (|24|) . The latter choice was made in the inset of 
Fig. 13. This yields a good description for N < 64. Two conclusions can be drawn from 
that: First, as expected from Eq. (|24j) , (X 2 (0)) carries the whole p-dependence so that ( 
is constant. This should be considered as a property of the bond-fluctuation model and 
not as a general feature because a bead-spring model can exhibit a different behavior ||61|| p|. 



Second, although the diffusion coefficients for N < 64 do not obey Eq. ( pq) [see Fig. |Tj] , the 
Rouse relation ( |24]) effectively remains valid. This is possible because the simulation results 
for t p can be written as r p = /(iV)[sin(p7r/2iV)] _2 and the dependence of the function / 
on N is just given by (iV(iV - l)D / B%,)~ 1 . Theoretically, N(N - l)D/R 2 e should represent 
the reorientation time of a monomer, (,b 2 /k B T, which is independent of N. In practice, 
iV(iV - l)D/R 2 varies from about 6 x 10~ 4 for iV = 16 to about 3.5 x 10" 4 for N = 64, 
i.e., by about 40%. For larger chain length, however, the residual iV- dependence of r p can 
no longer be lumped in a prefactor, although iV(iV — 1)D/R 2 decreases to about 7 x 10~ 5 



2 This point was not stated explicitly in Ref. |6T]]. However, since roughly (X 2 (0)) ~ p~ 2 ' 2 and 
t p ~ P~ 2 i n this case, Eq. (|24] ) would imply C = Cp ~ p ' 2 . Similar observations were also 
made in neutron scattering experiments on polyisobutylene . A quantitative description of the 
experimental results was only possible if both local stiffness and an increase of £ with increasing 
mode index p were taken into account. 
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for iV = 512 so that it is 1 order of magnitude smaller than for iV = 16. Nonetheless, r v 
depends more strongly on iV than NiN — 1)D/B% for small p and does not seem to be 
any longer a function of p/N alone. Even if one changes the prefactor to shift the data for 
N = 128 onto those for iV = 512, one can guess from Fig. [H| that such a collapse is only 
possible for p/N > 0.03, whereas t p for N = 512 becomes larger than that for N = 128 
when p/N < 0.03. Since scaling with p/N was possible for (JT 2 (0)), which probes the static 
properties of a chain only, Fig. suggests that the friction coefficient ( becomes p- and 
iV-dependent. 

Finally, Fig. H shows the variation of r ce [defined by Eq. ([H])] and of D with chain 
length. The ordinates were scaled by the chain length dependence expected from the Rouse 
model, i.e., r ee /N 2 and ND [see Eqs. (^2|) and (P5f)]. The figure illustrates that neither r ee 
nor D behave Rouse-like for small N or reptation-like for large N. A possible explanation for 
this finding is given by a recent reanalysis [|36| of the crossover scaling from the semi-dilute- 
solution to the melt state of the bond- fluctuation model fl3"3| . This new analysis focused on 
the diffusion coefficient, but substantially varied chain stiffness, which had not been done 
before. The important outcome for the present work is that the diffusion data of Fig. |H just 
lie in the crossover regime between Rouse dynamics, which is realized by flexible chains in 
dilute solution, and reptation dynamics which is most pronounced for stiff chains in the melt. 
Apparently, chain stiffness reduces the tube diameter and thereby amplifies the reptation 
behavior. The same conclusion has also been suggested by a recent molecular-dynamics 



simulation of the Kremer-Grest model ||54|| . Since the diffusion data belong to the crossover 
regime, this should also hold for r ee in Fig. [14| so that neither r ee ~ iV 2 nor r ee ~ iV 3 can 
be observed. If one just fitted the last two chain lengths (N = 128, 512) by a power law, 
one would obtain r ee ~ iV 3,23 . This is similar to experiments where one often finds that the 
relaxation time scales as iV 3 - 4 [§||. Of course, the present simulation data cannot prove, 
but only provide an indication that the bond-fluctuation model approaches (the possible 
asymptotic behavior) r ee ~ iV 3 from below with an exponent larger than 3, as suggested 
theoretically @,0,^1|]. These theoretical studies show that the crossover to the reptation 
asymptote is postponed to very large N. To observe this asymptote one has to study either 
much longer or (eventually) stiffer chains which seem to show a clearer signature of reptation 



for accessible chain lengths [136,54 



Nevertheless, one can use the simulation results for r ee to estimate the entanglement 
length N e by requiring that r e corresponds to the relaxation time of a chain with end-to-end 
vector d\ = R 2 (N e ). To read off N e , we fitted r ee in the range of r e (= 92114), i.e., for 
16 < iV < 64, by a power law. This gives r ee = 24.71iV 2 - 31 so that N e ~ 35. This value 
is compatible with our previous estimate from d\ = Rl(N e ), i.e., N e 38 (see Sect. |111 A| ), 
so that on average N e ~ 37. This result agree very well with that obtained by the above 
mentioned reanalysis of the crossover scaling [SH] . 



IV. CONCLUSIONS 



The present paper summarizes an attempt to extend previous work [|12|,[34],[35[ on the 
dynamics of long chain polymer melts towards larger chain lengths and longer times. The 
aim was to sufficiently enter the reptation regime to perform a careful test of the theory. 
However, it turned out that the crossover from non-entangled to entangled dynamics is 
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very gradual. Even the longest chain, N = 512, which is about 13 times larger than the 
entanglement length N e , can only be considered as slightly entangled when analyzed in terms 
of the asymptotic power laws of reptation theory. This finding is expected from a theoretical 



analysis of the Evans-Edwards model [|30|-p3|] and molecular dynamics simulations j56| , |57 



For N = 512 we find evidence that there is an entanglement time r e , above which the 
mean-square displacements of inner monomers, gi{t) [~ g§ (t)] and g2{t), and that of the 
center of mass, ^(t), increase as ~ t 1 / 4 and ~ t 1//2 , respectively. From these quantities the 
tube diameter d T and r e can be estimated: g?t ~ 18.6 and r c « 92114. Using these values 
and requiring Rl(N e ) = d\ and T ee (N e ) = r c , the entanglement length can be read off of the 
simulation data: N e ~ 37. The quoted values for <i T and r e are averages derived from g 2 and 
(76- The results from g 2 and g§ were rather close, whereas those from g% sometimes deviate 
substantially and thus appear less reliable (perhaps due to preasymptotic effects, as pointed 
out in ||33|| ). Of course, there are still uncertainties connected with these numbers, especially 
with g?t, since the prefactors are unknown and N = 512 is not really asymptotic. Here, the 
important point is that they are consistent with one another and with other independent 
studies. For instance, the value N e ~ 37 agrees very well with that found in [ ft6| . 

In addition to mean-square displacements evidence for reptation-like dynamics can be 
found in orientational correlation functions. Bond vectors in the innermost part of a chain 
decorrelate from their initial orientations more slowly than bonds close to the ends. This 
difference is present for all chain lengths, but becomes much more pronounced if N ^> iV c 
and t c <^ t ^ r ee . If there was a confining tube around a chain, such a behavior would 
be expected because reorientations of inner monomers are strongly hindered as long as 
the geometric constraints exist (i.e., for t < r ec ), whereas chain ends always take part in 
tube renewal and thus dissolve the confinement immediately. The simulations suggest that 
sensible indicators of this behavior are a late-time power law decay of the bond-vector 
correlation function and a modulated (eventually two-step) relaxation of the end-to-end 
vector or Rouse modes correlation functions as well as the ratio g^it) / gi(t) , as proposed 
in |]31|-p3|. The challenging problem consists in determining these functions with sufficient 
statistical accuracy in the relevant range of chain lengths and times. 

For short chains (N < N e ) the Rouse model represents a viable, though not perfect 
approach. Deviations are observed particularly on short length scales along the polymer 
backbone. This is to be expected because the microscopic potentials of the polymer should 
strongly influence the local dynamics. In general, these potentials cannot be represented 
by a flexible concatenation of Rousian harmonic springs, but introduce some local stiffness. 
If stiffness is taken into account, one can obtain a (almost) quantitative description of the 
static correlation (Xp(Q)) and of the p-dependence of the Rouse modes' correlation time 
t p , which scales as r v ~ (Xp(0)). Whereas the former property could be a general feature, 
the scaling r v ~ (X 2 (0)) is certainly a property of our model. Simulations of other models 



|61~| and neutron scattering experiments on polyisobutylene |52| suggest that the monomeric 



friction coefficient can be mode dependent. Such findings are plausible because one would 
in general not expect that static properties can uniquely determine the dynamic behavior. 

Another characteristic dynamic property of the Rouse modes is the progressive non- 
exponential character of the relaxation with increasing mode index p. This behavior is not a 
special feature of the bond-fluctuation model, but is also observed for a bead-spring model 
|B"T| and in chemical realistic simulations of polyethylene [ . These deviations from a pure 
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exponential decay are important if one wants to calculate the relaxation of quantities, for 
which large modes and short times also contribute significantly. An example is provided by 
the bond- vector correlation function (see Fig. [|). If one inserts the simulated Rouse modes 
instead of Eq. (|j"6|) into Eq. fll5D , a perfect description of the correlation function is obtained. 
Since this description only requires the orthogonality of the Rouse modes, it shows that the 
melt dynamics does not introduce (significant) correlations between different Rouse modes, 
i.e., (X p (t) ■ X q (0)) oc 5 pq . We believe that this property is general, since it is also observed 
for other models f61| . 

A further important feature is that there is no significant regime of chain lengths where 



the Rouse prediction D ~ 1/N holds (see inset of Fig. [nj). Furthermore, the related 
prediction gs(t) ~ t (for t > 1 in our units) is also not valid, not even for N = 16 (see 
Fig. H). If g 3 (t) ~ t were true, the log- log plot g^/t 3 ^ versus t should exhibit a straight line 
with slope 1/4. However, we find a curved behavior spread over several decades. Our results 
(Fig. 13) imply that there are about two decades in chain length (10 < iV < 10 3 ) which are 
neither fully described by the simple Rouse model nor by asymptotic reptation theory. A 
more complete theoretical description is thus called for. 
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TABLES 



N 


Rl 


R l 


(Ai) : <A 2 ) : (A 3 ) 




7"ee 


Tl 


16 


136 


23 


12.42 : 2.78 : 1 


391 


15180 


14519 


32 


293 


49 


12.30 : 2.74 : 1 


454 


78509 


76879 


64 


607 


102 


12.17 : 2.71 : 1 


430 


367863 


384993 


128 


1314 


217 


11.88 : 2.67 : 1 


445 


2387464 


2558552 


512 


5348 


885 


11.87 : 2.65 : 1 


461 


214679565 


233045839 



TABLE I. Chain length dependence of the size of a chain and of its relaxation times: R% = 
mean-square end-to-end vector, R^ = mean-square radius of gyration, (X a ) = mean eigenvalue of 
the gyration tensor in spatial direction a (= 1, 2, 3) [see Eq. ([]])], Tb = relaxation time of the bond 
vector defined by 4>b(T~b) = 1/e, r ee = relaxation time of the end-to-end vector defined by Eq. (21), 
ti = relaxation time of the first Rouse mode defined by (pb{Tb) = 1/e. 
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FIGURES 



FIG. 1. Mean-square end-to-end vector, R 2 (upper symbols), and radius of gyration, R 2 (lower 
symbols), plotted against chain length N. Asymptotically, both quantities are expected to behave 
as R 2 = Cool 2 (N - 1) and R\ = Cool 2 (N - l)/6 for dense melts. Here, I = (I 2 ) 1 / 2 (= 2.64) is the 
mean bond length and Cqq the characteristic ratio at infinite chain length. The value of Goo — 1-56 
was obtained by an extrapolation of R 2 /(N — l)l 2 to N — > oo using N > 100. Both R 2 and R 2 
enter the asymptotic regime if N > 100 (for R 2 this is hard to see on the scale of this figure; 
see Fig. ||). Furthermore, the figure shows that the present data are consistent with previously 
obtained results (on smaller systems) pi, 35]. 



FIG. 2. Eigenvalues, {X a )/(N — 1) (a = 1, 2, 3), of the gyration tensor Q [Eq. (|])] plotted versus 
chain length N. If the shape of a polymer were spherical, one would have (Ai) = (A2) = (A3). The 
result (Ai) S> (A2) > (A3) shows that it rather resembles a flattened ellipsoid ("prolate" chains). 
The ratio of the eigenvalues ((Ai) : (A2) : (A3) = 11.87 : 2.65 : 1 for N = 512) is fairly close to that 
expected for random walks (12.03 : 2.72 : 1 [45]). As anticipated from Eq. (|2|), the eigenvalues sum 
up to the radius of gyration calculated independently by R 2 = Yln=i(Rn ~ ^cm) 2 /^i where R n 
and .Rem are the position vectors of the nth monomer and of the center of mass, respectively. 



FIG. 3. Rescaled static autocorrelation of the Rouse modes, 8N(N — l)(X 2 (0)) /R 2 , versus 
p/N for all chain lengths studied (R e = end-to-end distance). This rescaling is motivated by 
Eq. (|4|). The dashed and the solid lines shows the theoretical predictions for fully flexible chains 
[Eq. (||)] and for freely rotating chains [Eq. (||)] where local stiffness due to the bond angle is taken 
into account. If p/N < 0.05 (corresponding to subchains larger than about 20 monomers), both 
predictions coincide and describe the simulation data well. In this case, Eqs. (|j) and (||) are well 
approximated by the power law p~ 2 (indicated by a solid line in the figure). However, if p/N > 0.3 
(corresponding to subchains smaller than about 3 monomers), the fully flexible (standard) Rouse 
model overestimates the correlation, whereas the (approximate) consideration of chain stiffness still 
provides a reasonable description. A magnification of the comparison between Eqs. (Q) and (|5|) is 
shown in the inset for N = 512. 

FIG. 4. Test of the orthogonality of the Rouse modes at t = [see Eq. (|j)]. Figure (a) shows 
(Xi(0) • X q (0)) (q = l,..., 127) for N = 128 (> N e » 37; N e = entanglement length). The inset 
magnifies the results for q = 2, . . . , 127. Cross-correlations between the first and the q = 2, 3 modes 
are about 1 to 2 orders of magnitude smaller than (X 2 (0)), but much stronger than those for q > 3. 
If N < N e , the behavior is qualitatively the same as shown in the inset. However, all correlations 
are more than 2 orders of magnitude smaller than the self-correlation for non-entangled chains. 
Figure (b) compares (X p (0) ■ X q (0)) of the largest mode for N = 128, (A m (0) • X q (0)}, with 
that for N = 16 (< N c ), (Xi^(0) ■ X q (0)). If N > N e , cross-correlations with the smallest Rouse 
modes are as large as the self-correlation, whereas they are at least an order of magnitude smaller 
if N < N e . 
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FIG. 5. Time dependence of gety/t^* for different chain lengths N. ge(t) is the cubic invariant 
defined by Eq. @. This function is predicted to behave similarly to the mean-square displacement 
of an inner monomer gi(t), but to exhibit the characteristics of reptation more clearly than gi(t) 
[30-33]. In the interval 200 < t < 5000, g&(t) increases roughly as i 1//2 . Then, it crosses over to 



a weaker time-dependence which is compatible with g§ (t) ~ i 1 / 4 for 4 x 10 5 < t < 2 x 10 6 . The 
intersection between these two asymptotic power-laws provides an estimate for the entanglement 
time T e « 100175 (solid vertical line). Using Eq. flT0|), g^T^/T^ 4 ~ d\/2>rl^ « 7, the tube 
diameter is approximately given by d-r ~ 19.3. Furthermore, the two dashed vertical lines indicate 
the values of the relaxation time, r ee , of the end-to-end vector, defined by Eq. (21), for N = 16 
and N = 512. 

FIG. 6. Time dependence of 92(i)/i 1 / 4 for different chain lengths N. gzit) is the displacement 
of an inner monomer relative to the chain's center of mass [Eq. @]. For times shorter than the 
Rouse, tr, or disentanglement time, Td, one expects gi (t) ~ g\ (t), whereas 92 (i) — * Rg otherwise. 
Due to Eq. ( |22] ) these times can be represented by the relaxation time, T ee , of the end-to-end 



vector [defined by Eq. (21)]. The latter behavior is borne out for t > r ee (vertical dashed lines 
shown for N = 16 and 512 only). As in Fig. |5| for g$ (~ g±), there is an approximate t 1//2 -increase 
for 200 < t < 5000 which crosses over to a weaker time dependence. This weaker dependence is 
(perhaps) compatible with gi ~ t 1 / 4 for iV = 512 over about half a decade (5 x 10 5 < t < 10 6 ). 
Using this interval the intersection with the extrapolation of the i 1//2 -behavior yields estimates for 
the entanglement time, r e , and the tube diameter, dx: T e pa 84053 (solid vertical line), f« 18.6 



[via Eq. @]. 



FIG. 7. Ratio of the mean-square displacement of the end monomer, g^ (t), and the inner 
monomer, gi(t), versus logarithm of time for different chain lengths N [see Eq. (||)]. The dashed 
vertical lines indicate the relaxation time, r ee , of the end-to-end vector for iV = 16 and iV = 128 [r ee 
is defined by Eq. (|2l|); see Table Q. The vertical dotted shows the value of the maximum, 94/91 = 2, 
expected from the Rouse model. This maximum should occur for t < tr. For a reptating chain, 



a maximum is still expected to occur at t ~ tr, but to be much larger than 2 [31 33] . The onset 
of this behavior is visible for iV = 128. In the limit t — > 00, both displacements become diffusive. 
So, 94/91 — > 1. Since 94 ~ 91 in general, the calculation of the ratio 94/91 eliminates the dominant 
variation with time and highlights the differences between both displacements. Thereby, statistical 
errors also become much better visible. For N = 512, they are so pronounced that a quantitative 
analysis was not possible. 
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FIG. 8. Time dependence of gz{t)/t 3 ^ for different chain lengths N. gs(t) is the mean-square 
displacement of the chain's center of mass [Eq. (H)]. Since reptation theory predicts a crossover 
from 53 ~ t 1 ! 2 for r c < t < tr to 53 ~ i for t > tr, 53(i)/t 3//4 should exhibit a minimum at f ~ tr 
and increase as t" 1 / 4 and t 1 / 4 for t < tr and i > tr, respectively. Qualitatively, such a behavior 
is observed for N = 512 only (solid lines with slope 1/2 and 1). The intersection point of both 
power laws yields an estimate for the Rouse time tr ~ 4.34 x 10 6 , whence di ~ 13 by virtue of 
93 ( t r) = d^/3. The minimum is preceeded by another power-law increase (empirically, 53 ~ 
if i < t c . For smaller chain lengths the power law is not as well pronounced. The intersection 
point of the power laws 53 ~ i ' 84 and 53 ~ i 1 / 2 yields another estimate of the entanglement time, 
r e ~ 268718, and of the tube diamater, ~ 20.3 [see Eq. (|l3|)]. 53 becomes diffusive in the 
range t ~ r ee [indicated by dashed vertical lines for iV = 16 and 512; r ee = relaxation time of the 
end-to-end vector, see Eq. (|2~i"|)]. 



FIG. 9. Scaling plot of the autocorrelation function, 4>h{t), of the bond vectors for different 
chain lengths versus t/r^. The relaxation time, Tb, was defined by (see Table |): c/>b(^b) = V e 
(dashed horizontal lines in the main figure and the inset). This scaling collapses the correlators 
onto a master curve if 0.25 < 4>b(t) < 1. The final decay, 4>h{t) < 0.25, depends on chain length. 
For N = 16, it is well described by the Rouse model, i.e., by Eq. ( |i~5|) (dashed curve in the main 
figure), whereas Eq. ( |i~5| ) decays faster than the simulation data for N > 32. Especially, the 
largest chain lengths (N = 128,512) exhibit a slow power-law decay (emprically, 4>b{t) ~ t -0 ' 35 ), 
as illustrated in the inset. The inset also indicates the Rouse behavior, </>b(i) ~ t~ a5 , expected for 
t < r c (r c = 92114 = average entanglement time of Figs. || and ||). 



FIG. 10. Survival probability (17), 4> n (t), of an initial bond vector l n (0) for N = 16 [panel 
(a)] and N = 128 [panel (b)]. (f) n (t) is plotted as a function of the bond number n. For example, 
n = 1 corresponds to the bond vector ii(0) connecting the first and the second monomer. Time 
is measured in units of the relaxation time of the end-to-end vector r ee (see Table |). The solid 
lines represent the prediction of Eq. (|l8|). To improve the statistics the simulation results were 
averaged over both halves of the chain - the chain is symmetric with respect to its ends. Therefore, 
results from n = 1 to n = N/2 are shown only. This was sufficient for the short chain, but not for 
N = 128. In this case, the bonds n and n + 1 were lumped additionally and the resulting average 
was plotted at the (non-existing) bond number n + 1/2. 



FIG. 11. Autocorrelation function of the end-end vector, (j} e {t), versus scaled time t/r ee for all 
chain lengths studied. The scaling time, T ee , is defined by Eq. (||): M T ee) = 0.298221 (dashed 
horizontal line). The vertical dashed line indicate the value of the entanglement time r e (= 92114 = 
average r ee from Figs. |5| and |6|) and of the Rouse time tr(A^ = 512) [~ 4.34 x 10 6 ] determined 
from Fig. ||. Both times are divided by T ee (N = 512) (see Table |). The solid line depicts the 
result of Eq. (f20|). The inset shows 1 — <pe{t) versus t/r cc for iV = 512 together with the power laws 
1 — 4> c (t) ~ t 9 / 32 and 1 — (pe(t) ~ t 3 / 8 expected from polymer mode-coupling theory for r c ^ t <ti tr 
and tr <C t <C trr (trr = renormalized Rouse time) [24). 
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FIG. 12. Correlation function of the Rouse modes, <j> p (t), versus rescaled time t/r p for N = 16 
[panel (a)] et N = 512 [panel (b)]. The scaling time, t p , is defined by: (f> p (r p ) = 1/e (dashed 
horizontal lines). For both chain lengths the relaxation of four representative mode indices p 
(= 1, ... ,N — 1) is shown. In addition, an exponential function (solid lines) is depicted. This is 
the behavior expected from the Rouse model [see Eq. fll6|)1. The vertical dashed line in panel (b) 
indicates the Rouse time tr(N = 512) [~ 4.34 x 10 6 ] determined from Fig. ||. tr is divided by the 
relaxation time, ri, of the first Rouse mode (see Table Q). 



FIG. 13. Variation of the Rouse mode relaxation time with the mode index p for different chain 
lengths iV (entanglement length N e ~ 37). The relaxation time, r p , is defined by the time value at 
which the correlation function of the Rouse modes has decayed to 1/e, i.e., 4> p {t p ) = 1/e. Motivated 
by Rouse theory the abscissa and ordinate are scaled according to Eq. (|23|) where Q? jk-^T is 
identified with R 2 JN(N-1)D (D = diffusion coefficient, see Eq. @; R 2 = mean-square end-to-end 
vector). This scaling leads to a reasonable collapse if N < 64, whereas the p-dependence for 
N = 128, 512 (^> N c ) is s-shaped. t p is larger than the Rouse prediction (solid line) if p/N is small 
(especially for iV = 512), but smaller than it for p/N — > 1. As for the static correlation (X 2 (0)) 
(see Fig. ||), the Rouse model provides a good description of the p-dependence if p/N < 0.05 (and 
N < 64), but overestimates the relaxation time for p/N > 0.3. The inset therefore tests whether 
t p can be written 2(X 2 (0))/3£> [see Eq. @] with {X 2 (0)) taken from the simulation (see 

Fig. |3|). This is possible for N < 64, but not for large (entangled) chains, such as N = 512. 



FIG. 14. Chain length dependence of the end-to-end vector relaxation time, r ee (main figure), 
and of the diffusion coefficient, D, of a chain (inset). r ee and D are determined from Eq. ( ^l|) and 
from the long-time limit of g% (= limj-^oo gs(t) /6t), respectively. Both quantities are scaled by the 
TV-dependences of Eqs. (|22|) and (|25| ) to highlight deviations from Rouse behavior. The diffusion 
coefficient is compared with results from previous studies [[34.35]. The solid lines in the main figure 
and the inset indicate the behavior, r ee ~ A^ 3 and D ~ iV~ 2 , suggested by reptation theory. The 
triangles (A) are the predictions of the Rouse formula r oe = {Nb) 2 C i /?,iT 2 k-QT = R 2 /3ir 2 D, where 
the simulation results are used for D and R 2 (= end-to-end distance). 
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